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Abstract 

Consequences resulting from the DD excess derived indirectly by the NA50 experiment in 
S+U and Pb+Pb collisions at y/s=19, 17 GeV, relevant for the identification of the QCD 
phase transition in these collisions, are discussed. The dependence of open and closed 
charm yields in Pb+Pb collisions on the number of participating nucleons (N) indicates 
non thermal charm production and J/^P dissociation, stronger than the absorption seen 
in any other elementary hadron. The J/'^ in central Pb+Pb collisions could originate 
dominantly from cc pair coalescence out of a hadronizing quark and gluon environment. 
Furthermore, the J/\Ef appears to be suppressed in S+U collisions at y/s=19 GeV, as 
opposed to current interpretations. A significant change in the {J/'^)/DD ratio as well 
as in the number density of kaons is observed above energy density e ~ 1 GeV/fm^, 
suggesting a change of phase at this energy density, and underlining the importance of 
direct open charm measurements. 



1. Introduction 



Quantum chromodynamics (QCD) on the lattice predicts a phase transition of confined 
hadronic matter into deconfined quark and gluon matter (called the Quark Gluon Plasma 
state - QGP) at a critical temperature ~ 150 MeV, respectively at energy density ec ~ 
1 GeV/fm^ [|l|. The order of the transition is parameter dependent 0. An investigation 
of relevant observables in heavy ion collisions |^ as a function of energy and/or the impact 
parameter^ of the collision could reveal this transition, through a discontinuous behaviour 
of many QGP signatures at the transition point. 

An anomalous suppression of the J / \& meson predicted to be a signature of Quark Gluon 
Plasma formation |^ has been measured to occur in the ratio of the J/\E' over the Drell Yan 
(DY) process in Pb+Pb collisions at ^/s=VJ GeV investigated as a function of transverse 



^Work supported by the Swiss National Science Foundation. 

^Provided that there is a unique assignement between impact parameter and QGP phase transition. 
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energy (Et) 0. In S+U collisions at v^=19 GeV and in the most peripheral Pb+Pb col- 
lisions, the {J 1'^)/ DY ratio agrees with expectations 0. The ratio {J/'^)/ DY is relevant 
for the investigation of the J/\l/ suppression, under the assumption that J/^ production 
in these collisions is a hard process. 

Recent measurements of the dimuon invariant mass m(/i"'"/i~) spectrum between the (f) 
and the J/\l/ mass (Intermediate Mass Region=IMR) revealed a dimuon enhancement 
above expectation, which is increasing with the number of nucleons participating in the 
collision (N) ||^. This enhancement can be understood as due to an excess of DD pro- 
duction| as suggested by several features of the data, e.g. the shape of the mass, rapidity, 
angular and transverse momentum distributions of the dimuons 0. The interpretation 
of the IMR enhancement as due to open charm is not unique though, because the open 
charm was extracted through a fit to the dimuon continuum and the DD meson signal 
was not directly identified. Recent work suggesting that the seen enhancement could be 
due to rescattering of D mesons in nuclear matter is not supported by the data [^]. 
An other possible interpretation is that the IMR excess could be due to thermal dimuons 
[y]. Non-perturbative effects are known to play a role in heavy flavour hadron produc- 
tion in elementary reactions showing up in deviations of data from perturbative QCD 
calculations |llO| , [TT| . Based on this fact, one could expect that theoretical investigation 
of non-perturbative effects and different reaction dynamics as in the plasma phase may 
result in an enhancement of open charm production in nuclear reactions over perturbative 
QCD expectations. 

If the total charm produced in heavy ion collisions indeed deviates from the perturbative 
QCD expectations for a hard process as suggested by the NA50 data, it follows that the 
{J/'^)/DY ratio is not the proper quantity for the search for the J/\l/ suppression as 
signature of Quark Gluon Plasma formation in nuclear reactions. It is only the ratio 
(J/\l/)/(total cc) that matters. We therefore investigate here first the dependence of the 
J/^ and the DD yields per collision on N. We further investigate the dependence of the 
{J/'^)/DD ratio on N and on the length of the nuclear matter traversed by the J/\l/, as 
well as the dependence of both charm and strangeness production on the initial energy 
density reached in the collision. 

2. N dependence of open and closed charm yields in Pb-|-Pb collisions at 
v^=17 GeV 

2.1. N dependence of the Drell Yan yield 
Calculation details 

The dependence of the Drell Yan (DY) yield per nucleus-nucleus collision in arbitrary 
units produced in Pb+Pb collisions at y/s=VJ GeV per N+N collision, on the transverse 
energy of the collision has been measured by the NA50 collaboration (figure 7 in ll^). In 



this figure it is shown that the theoretically expected DY yield, assuming DY production is 



■^With we denote the number of D and D hadrons which were simuhaneously found within the 

acceptance of the NA50 experiment ( DD = {D + D)acc/'^ )■ 
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a hard process, does partly reproduce the measured one; the deviations at low transverse 
energy are understood to result from properties of the lead (Pb) nucleus, in particular 
from the different radii of the proton and neutron distributions in Pb |T^. It is therefore 



justified to use the theoretically calculated Et dependence of DY yield per coUision, to 
estimate the dependence of J/\E' and DD yields on the number of participating nucleons 
in the collision. 

However, since the deviations of the very low statistics measurement of DY yield at high 



Et, from the theoretical estimated DY yield seen in figure 7 of |T2| cannot be understood, 
it would be important to measure the DY yield per collision with high statistics in the 
high Et region. In this way the last drop in the {J/'^)/DY ratio above N ~ 360 can 
be experimentally verified. 



2.2. N dependence of the DD yield 



Calculation details 

The NA50 collaboration observed an excess (E) of the measured over the expected DD/ DY 
ratio in S+U and Pb+Pb collisions at -^5=19, 17 GeV, which increases with the number 
of participants N (figure 12 and table 4 in [0). If we fit the S+U and Pb+Pb E points of 
the above figure to a function / = c ■ iV", we find that the excess is increasing with N as 
jY(a=o.45±o.ii) (^^Degrees Of Freedom=1.7, D0F=7). The dependence of the excess 
E of the DD/DY production in S+U collisions at y^=19 GeV and Pb+Pb collisions at 
y^=17 GeV over expectations, reflects the dependence of the DD/DY ratio. This 
results from the fact, that all other quantities involved in the definition of E P,^], do not 
depend on A^. Therefore the A^ dependence of the DD production yield is given by the 
A^ dependence of the quantity 

njjjj = E* udy ~ (DD/DY) * udy (1) 

where n^yj^, udy denote the yields of DD and DY per collision in arbitrary units. The 
arbitrary units are due to the fact that NA50 did not published absolute yields per colli- 
sion of the DY and DD separately, corrected for losses due to e.g. acceptance, as a 
function of N, Et- We suggest that it would be important to do so. 

The DY yield used for the above calculation has been extracted from the theoretical curve 



shown in figure 7 in [12|, at the transverse energy (Et) points in which the DD excess 
factor E has been measured. The Et points corresponding to the excess factor E were 
extracted, by interpolating between the Et values given in table 1 of [|l^ as a function 



of the mean impact parameter, at the values of mean impact parameter for which the 
factor E has been measured (listed in table 1 of 0). For the most central and the most 
peripheral points, for which no mean b are given in table 1 of [^j , we estimated the values 
of b, from the values of A^ as a function of b for Pb+Pb collisions calculated by [|1^,|15[. 
These calculations fl^ agree with the values (N,b) estimated by NA50, when compared 
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in their common range. 

Though the DD measured by NA50 represents the joint probabihty that a D and a D are 
both found in the experimental acceptance, the N dependence of it, is expected to be the 
same or similar to the N dependence of the sum D + D, respectively of the D, and of the 
D hadrons. It is however possible that the above N-dependences deviate from eachother, 



in case of a higher than 1 charm pair multiplicity per event with charm (see [0 for a 
numerical estimation of the charm multiplicity per event in central Pb+Pb collisions at 
158 A GeV). Therefore, for clarity, the N dependence of the total extrapolated D + D 
yield in nuclear collisions should be estimated including acceptance corrections by NA50. 

Results and discussion 

The resulting DD yield in arbitrary units (figure |lD increases as ^^("=1-70=1=0.12) (^y^ j jjQp 
= 2.5, D0F=7). This N behaviour indicates that DD production in Pb+Pb collisions 
at y/s=17 GeV, did not establish yet equilibrium, in which case a proportionality with 
-assuming N to be proportional to the volume of the source is expected {a= 1). 
This appears justified because the temperature in the collision zone -assuming local ther- 
malization of light particles- drops with time and the mean temperature expected to be 
reached in these collisions of the order ~ 10^ MeV, is much lower than the mass of charm 
quarks and/or charmed hadrons. 



2.3. N dependence of the J/vp yield 



Calculation details 

In the following we estimate the J/^P yield per collision as a function of A^, at the same 
A^ values where the DD was measured. The A^ dependence of the J/\l/ yield per collision 
is given by the N dependence of the quantity : 

nj/^ = {{J/^)/DY)*nnY (2) 

where nj/ij,, udy denote the yields of J/\l/ and DY in arbitrary common units. The 
{J/'^)/DY values have been extracted from figure 4 of |^ at the values where the E 
factor has been measured, interpolating between the different points. 



Results and discussion 

The resulting J/\l/ yield per collision produced in Pb+Pb collisions in arbitrary units (fig- 
ure|) increases like iV("=o-70±o.o4) {^^^^lDOF=lA'i, D0F=7). This N dependence indicates 
an increasing J dissociation with higher centrality. The strength of this dissociation 
as measured by the a parameter, is higher for the J/\& as compared to any other hadronQ 

^ The assumption N V, is based on the observation that the freeze out volume V of the particle source 
is found to be proportional to N p7[ |. 

^Deuterons have an even smaller a parameter, but they are not elementary hadrons and are weekly bound 
(see discussion in [TtI]). 
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produced in these collisions for example as compared to antiprotons. For the latter, a 
large annihilation cross section with baryons is expected and there is indeed experimental 
evidence that they are absorbed with increasing centrality in Pb+Pb collisions (q;(p)=0.80 
± 0.04, (xVi5OF=1.0, D0F=3) at y=3.7, pt=0 fT^Q). 



The J/^ multiplicity as a function of N extracted with an other method [|T8[, agrees with 
the here presented results within the errors. 



3. The {J/'^)/DD ratio in nuclear collisions 

3.1. The N dependence of the {J/'^)/DD ratio in nuclear collisions 
Calculation details 

Assuming that the IMR excess is due partly or solely to open charm, allows us to search 
for an anomalous suppression of J/\E' as compared to the open charm production. The 
N dependence of the {J /'^) j DD ratio in Pb+Pb and S+U collisions at i/s of 17 and 19 
GeV, estimated as: 

{J/^)/DD ~ {{J/'^)/DY)/{DD/DY) ~ {{J/'^)/DY)/E (3) 

in arbitrary units due to the E factor in equation (3), is a decreasing function of (fig- 
ure H). The {J/'^)/DY in S+U collisions was taken from [|19| , ^] . Note that possible 



deviations of the DY yield from its theoretical calculation (as seen in figure 7 in |T^), 
do not drop out in the {{J / / DY) / {D D / DY) ratio shown here, because the DD/DY 
-unlike the {J /'^) / DY- was calculated by NA50 not using the minimum bias theoretical 
DY yield values but the measured ones. 

In order to show the influence of the very last drop of {J/'^)/DY, on the {J/'^)/DD 
ratio, the {J/'^)/DY ratio divided by A^o-45±o.ii jg plotted function of N in figure H. 
This quantity resembles the {J/'^)/DD ratio 

( J/^) /DD ~ (( J/^) /DY) /ArO-45±o.ii 

in arbitrary units, because is the found N dependence of the DD/DY ratio. 

The open points of figure |^ are extracted from the 'minimum bias analysis' results of 
figure 4 in §[]. The closed points show the {{J /-^)/ DY)/N^-^^ calculated here, at the N 
values where the DD excess factor was measured. 

Results and discussion 



^We extracted here the a parameter for p, after quadratically adding the statistical and a 5% systematic 



error 

7) 



Minimum bias analysis' in NA50, means that the DY for the (J/'^)/ DY ratio, was determined using 
the theoretically estimated DY yield per collision as a function of Et and the measured dN/dEx vs Et 
spectrum of minimum bias trigger events (see |l2| ). 
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The {J/'^)/DD ratio in Pb+Pb collisions found as shown in equation (3), decreases 
with N as A^(°=-o.79±o.i4) {^^^ijjQp ^ 3 3^ D0F=7), respectively like N^»=-^-^7±o.u) 

{X^/D0F=1M, D0F=6) when the first point is not fitted^. The {J/^)/DD ratio in 
S+U collisions (figure D decreases with as A^("=-o-62±o.22) ^^yDOF=0.Q9, D0F=3). 
The {J/'^)/DD found as shown in equation (4) when fitted to the function A^" until 
N=380, gives A^{"=-i-07±o.o7)^ (x/£'OF=0.93, D0F=7). 

If the J/^ is completely dissociated in a quark gluon plasma and is formed later mainly 
through c and c quark coalescence, we expect that the N dependence of the ratio {J j^)/ DD 
-rather than the {J/'^)/{DD)'^- reflects the N dependence of the volume of the charm 
environment |10[- This is due to the expectation that because of the very low cross sec- 
tion of charm production at these energies, there is most often just one cc pair per event 
containing charm, whatever N. Then the probability to form a J/\E' from coalescence is 
proportional to {J /'^)/DD and inversely proportional to the volume of the particle source 
-made up by uuddss quarks and gluons- within which the c and c quarks scatter. Assum- 
ing this volume is proportional to N (see footnote 3), one would expect that {J/'^)/DD 
decreases as A^~^, as actually observed. 

In this case, one can use the {J/^)/DD ratio to extract the absolute value of the volume 
of its environment with a coalescence model. The 'charm' coalescence volume would re- 
flect partly the QGP hot spot volume and partly the hadronic source volume from which 
hadrons with charm and anticharm can also form a J/\E'. If the absolute yields per colli- 
sion of J/"^ and DD as a function of N, needed for this calculation would be published 
by NA50, the charm coalescence volume could be calculated. 

Figure |^ suggests that the coalescence picture could hold for the full N range of Pb-|-Pb 
collisions up to N=380. Obviously, it would be better to use the DD data themselves 
instead of the N parametrization, if high enough statistics would be available. 

On the other hand, if the multiplicity of charm quarks is high enough that often more than 
1 charm quark pair per event with charm is produced, then it is the ratio {J/'$)/{DD)^ 
which is expected to be inversely proportional to the volume of the charm source (this is ex- 
actly the case if d coalescence out of p and n is investigated in a baryon rich source). The N 
dependence of the {J j'^)/ {DD^ ratio, which would be relevant in the above discussed case 
is n("=-2-26±o.^s), {x'^/DOF=2.5, D0F=7) respectively ]^(.'^=-3-i±o.24.) ^ {x^/D0F=1.2, 
D0F=6) if the first point is not fitted. The question on the absolute multiplicity of 
charm in nuclear reactions, should be answered by experiment. 

3.2. The L dependence of the {J /^!) / DD ratio in nuclear collisions 
Calculation details 



®The first point of the DD/DY enhancement factor E lyes significantly above the N" function fit to the 
E distribution (figure 12 in p). 
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The two distributions of {J /'^) / DD ratio for S+U and Pb+Pb collisions in figure are 
measured at different energies, and therefore they cannot be compared in terms of their 
absolute yields but only with respect to their shapes. In order to compare their absolute 
yields, the data from figure 5 of |^ will be used. There the {J/'^)/DY ratio in p+A, 
S+U and Pb+Pb collisions is shown as a function of L, all normalised to the same energy 
(-ys=19 GeV) and corrected for the isospin dependence of DY production. The param- 
eter L, is the length that the J /'^ traverses through nuclear matter. In order to convert 
figure 5 of to the {J /'^) / DD ratio as a function of the L parameter, the {J/'^)/DY 
ratio data points have been divided by the E factor as described in equation (3). The 
correlation of the L parameter with N and b for Pb+Pb collisions, has been estimated 
using the theoretical calculation of pi . 



Results and discussion 

The L dependence of the {J /'^)/ DD ratio in arbitrary units in p+A, S+U and Pb+Pb col- 
lisions calculated here, is shown in figure |^, together with the {J/'^)/ DY ratio published 
T9|j20[| . The closed points show the {J/'^)/DD ratio in S+U and Pb+Pb collisions 



m 



extracted as indicated in equation 3. The open squares and circles show the {J/'^)/DY 
ratio in S+U and Pb+Pb collisions from [|19| , p0| . The open stars, show both the L depen- 
dence of the {J/'^)/ DY as well as the L dependence of the {J /'^) / DD in p+A collisions 
which are the same, since the factor E has the value 1 for the latter. 



The J/\E' over the DD production investigated as a function of the volume through which 
the J /'^ traverses (V ~ L^) is suppressed as compared to the shape of the exponential fit 
going through the {J/'^)/ DD p+A data, in both the S+U and Pb+Pb collisions at all L 
points (figure H). 

The energy density of the lowest S+U point has been estimated to be ~ 1.1 GeV/fm^ 0, 
which is comparable to the predicted critical energy density for the QGP phase transition 
of ~ 1 GeV/fm^. A similar energy density of 1.2 GeV/fm^ has been estimated to be 
reached in the most peripheral Pb+Pb collisions measured by NA50. 



In the following we investigate the initial energy density, rather than only the volume of 
the particle source (V ~ L^), as a critical parameter for the appearance of the QGP phase 
transition. 



4. The e dependence of charm and strangeness in nuclear collisions 



4.1. Charm 



Calculation details 
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We estimate here the {J/'^)/ DD ratio as a function of the energy density e. For this pur- 
pose we use part of the data shown in figure 7 in p. There the ratio of {{J /'^) / DY) measured 
over the {{J /'^) / DY) expected is shown. The \{J l"^) / DY) expected, is taken to be the expo- 
nential fit seen in figure^, which represents the 'normal' J /'^ dissociation (i.e. understood 
without invoking QGP formation). Dividing these data points by jv'^'*^^'''^^, and normal- 
ising the distribution of S-l-U and Pb+Pb points to the p-|-p and p-|-A data as in figure 
^, we estimate the {{J/'^)/DD) ratio over the expectation expressed by the above men- 
tioned exponential curve, which fits the {{J j'^)/ DD) data points for p-|-p, p-|-d and p-|-A 
collisions. 

Results and discussion 

The result of this calculation is shown in figure ^ in logarithmic scale and in figure |^ (a) in 
linear scale. It demonstrates a deviation of the {J/'^)/DD ratio both in S+U and Pb+Pb 
collisions, from the p+p and p+A expectation curve, occuring above e ~ 1 GeV/fm^. The 
logarithmic scale is shown to reveal small changes in the slope of the {{J/'^)/ DD) distri- 
bution as a function of e, appearing at e ~ 2.2 and 3.2 GeV/fm^. 



4.2. Strangeness 

Figure |^ compares the two QGP signatures of J/\E' suppression and of strangeness en- 
hancement. For this purpose we represent all data points as a function of the estimated 
energy density. Note that the energy density as critical scale variable, has the advantage 
that unlike the temperature, it is defined irrespective of whether equilibrium is reached 
in the collisions studied. 

Figure ^ (b) shows the multiplicity of kaons per event {K^ , but also some data scaled 
to are shown) divided by the effective volume of the particle source at thermal freeze 
out in the center of mass frame, as a function of the initial energy density. The effective 
volume represents the part of the real source volume, within which pions are correlated 
with each other (the so called 'homogenity' volume in the literature |^). The effective 
volume is smaller than but proportional to the real source volume. For a more precise 
calculation of the freeze out source volume a detailled model is needed. Here we estimate 
the effective volume at thermal freeze out Vthermai based on measurements. The (smaller) 
effective volume at chemical freeze out Vchemicah is not experimentally measured, we give 
however an estimate of the ratio Vehemicai/ythermai- Note that we compare the kaon data 
without rescaling for the different energy between AGS and SPS[|. 

Calculation details 

The effective volume V of the particle source has been estimated in the center of mass 

^ The total multiplicity in p+p collisions, increases by a factor of '--^ 5, from 11.1 to 158 GeV per 
nucleon jH]. 
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frame, assuming cylindrical shape of the source: 

V = {-K- Radiuslyn^^^^) ■ LengthcyUnder ^ 

V = (n ■ A- i?Le) ■ iVT2 ■ -Rlong) 



where Rside is a measure of the transverse radius and Riong is a measure of the longitudinal 
radius of the particle source, and the factors 4 and \/l2 arise from the definition of Rside, 
The Rside and Riong values for central Au+Au collisions at 10.8 A GeV and for 



R, 



long 
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central Pb+Pb collisions at 158 A GeV have been taken from [^. We dont use the more 
elaborated estimation of the homogenity volume given in because the Roi component 
is not given in Based on the data of |^ we estimated the effective volume of the 
source at thermal freeze out in central Au+Au collisions at 10.8 A GeV {V ~ 1949/m^) 
and central Pb+Pb collisions at 158 A GeV {V ~ 6532/m^). The effective volume in- 
creases by a factor of 3.35 from AGS to SPS energy. 



Based on the temperature at thermal and chemical freeze out which has been estimated 



from measurements using thermal models [25|, and the above estimated volumes at ther- 



mal freeze out, we can further estimate the volume at chemical freeze out. For this we 
assume that the relation V ~ T~^, which holds in the universe for massless particles 



in thermal equilibrium and for adiabatic expansion p6|, holds approximately for heavy 



ion collisions at AGS and SPS energy. Then from the temperature values at thermal 
and chemical freeze out given in ||2^ averaged over all models, we find that the ratio 

Vchemical/Vthermali-^GS Si+Au 14.6 A GeV) =0.45 and Vchemical/Vthermal{SPS Pb+Pb 158 

A GeV) = 0.28. Using the volume at chemical freeze out as estimated above, would 
stretch the K/V ratio in figure ^ (b) between SPS and AGS, by a factor ~ 1.6 apart. 
We dont use these values in figure |^, because the above calculation is model dependent, 
e.g. the assumption of massless particles is not met, while the assumption of thermal 
equilibrium may not be true. 



The ratio K/N is expected to be proportional to the number density of kaons ~ K/V, 
(V=volume), assuming V ~ N (for justification of this assumption see footnote 3, ||23| 
and |2^). Based on this expectation, we estimated here the K/V ratios from the K/N 
ratios, by normalising the K/N ratios to the K/V value of the most central Au+Au events 
of E866 respectively Pb+Pb events of NA49, for which the value of the volume has been 
estimated above. 



The kaon data from Au+Au collisions at 11.1 A GeV (E866 and E802 experiments) ||28| 



and from Pb+Pb collisions at 158 A GeV (NA49 experiment) are kaon multiplicities 
extrapolated to full acceptance. Therefore NA49 and E866 data are absolutely normalised. 
We estimated K/N from the NA49 experiment using the kaon multiplicities from and 
the number of wounded nucleons from as available^, otherwise we used the N esti- 
mated from the experimental baryon distribution ||29|| . 



We take N equal to the number of wounded nucleons for NA49, because it is used in all other experi- 
ments presented here, and allows for a straightforward geometrical interpretation of N. 
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The data from NA52 and WA97 have been measured in a small phase space 
acceptance and have been scaled here arbitrarily, in order to match the NA49 data in fig- 
ure 0. This scaling is justified since all NA52, NA49 and WA97 measurements are kaons 
produced in Pb+Pb collisions at 158 A GeV, and 'extrapolates' the NA52 and WA97 
data to the NA49 full acceptance multiplicities, allowing for comparison of the shapes of 
the distributions. It is assumed, that the N and e dependence of kaons does not change 
significantly with the phase space acceptance. 

In order to calculate the energy density we have performed the following steps. The energy 



density for all colliding systems has been estimated using the Bjorken formula |33] and 



data given in p8| , |33| , |3^ . The transverse radius of the overlapping region of the colliding 
nuclei is found as: Rtrans = 1.13 ■ where N is the total number of participant 

nucleons. The formation time was taken 1 fm/c ||32| . 



For the E866 experiment, lacking Et values, (but with measured 'Ei forward), we used 
instead of {dEx / dvi)ycm, (|/cm= mi drapidity) the total energy of the nucleons participat- 
ing in the collision {Etot,part = Nprojectiie,participants " Ebeam) , assumiug the proportiouafity 
[dEx/ drj)ycm ~ Etot,part, and we further normalised the results in such a way that the max- 
imum energy density of our estimate matches the absolute value of the maximal achieved 
energy density in the most central Au+Au events at this energy of 1.3 GeV/fm^, given 
33[. 



m 



NA52 measures Et near midrapidity (y ~ 3.3). These values were used to estimate the 
energy density and the results have been normalised to the maximum energy density 
reached in Pb+Pb collisions at the same centrality of emax=3.2 GeV/fm^, extracted by 



NA49 in |24[. Parametrizing the dependence of the energy density on the number of par- 
ticipants found from the NA52 data as described above, we estimated the energy density 
corresponding to the N values of the WA97 and the NA49 kaon measurements, given in 



!]],B9^. Data from S+S collisions taken from 29 and 35,34 are also shown. 



To estimate the systematic error on the energy density found with the above methods, we 
calculated the energy density in Pb+Pb collisions at 158 A GeV, using the VENUS 4.12 
Q event generator. We estimated with VENUS the {dET/dr])ycm at ?/cm=midrapidity 



and the number of participant nucleons and used them to find the energy density from 
the Bjorken formula [Q. The deviation of the energy density calculated with VENUS 
{dET/drj)ycm from the energy density found using the NA52 transverse energy measure- 
ments is < 30% of the latter. The deviation of the energy density calculated with VENUS 
dEx/dr] from the energy density found using the total energy of the participant nucleons 
and of the newly produced particles, (which is similar to the method used to estimate the 
energy density for the AGS data), over the latter energy density, is at the same level. 



In this context, it appears important for a more precise comparison of data as a function 
of e, that experiments publish together with the number of participants also the dEx/dr] 
at midrapidity for each centrality region, for both nucleus+nucleus and for p+p collisions. 
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estimated by models or measured if available (e.g. in NA49). 



Results and discussion 



Figure |^ (b) suggests that kaons below e ~ 1 GeV/fm^ did not reach equilibrium, while 
this seems to be the case above. Indeed kaons produced in Au+Au collisions at 11.1 
A GeV and in very peripheral Pb+Pb collisions at 158 A GeV [ P7| , P5| , P^ , increase 
faster than linear with N, indicating non thermal kaon production, while they increase 
nearly proportional to N above e ~ 1 GeV/fm^ [^,^,^. The connection of strangeness 
equilibrium and the QGP phase transition has been discussed e.g. in There it is 

shown, that strangeness in heavy ion collisions is expected to reach equilibrium values if 
the system runs through a QGP phase, while this is less probable in a purely hadronic 
system. 



Figure |^ demonstrates that both the J/\l/ and kaon production exhibit a dramatic change 
above the energy density of ~ 1 GeV/fm^. While the equilibration of strange particles 
as suggested by their ~ A^^ dependence above 1 GeV/fm^, could in principle also be due 
to equilibrium reached in a hadronic environment, the combined appearance of this effect 
and of the {J/'^)/DD suppression at the same energy density value is a striking result, 
indicating a change of phase above ec=l GeV/fm^. 



The expectation for the shape of the J/\E' suppression as a function of energy density are 
three succesive drops of the J/\l/ P0|,^; a drop by ~ 8% ||TH] due to ip dissociation, a 



drop by ~ 32% |jT8[ due to the Xc dissociation and a drop by ~ 100% due to the J/\l/ 
dissociation. All these without taking into account regeneration of J/\E' through other 
processes. These can be e.g. coalescence of charm quarks or J/^ not travelling through 
the plasma. The ip feeds only 8% of the total J/\l/'s and can therefore hardly be observed 
as a break in the J/\l/ production. 



The absolute value of the energy density e and therefore of the N values at which these 
changes could be observed is not exactly given by the models. The critical energy densi- 
ties for the dissociation of the states \1/ , Xc and J/"^ could even be so near to each other 
that no clear multistep behaviour is seen in {J/'$)/DD. 

Figure suggests that the breaks in the {J/^!)/DD ratio at e ~ 2.2 and 3.2 GeV/fm^ 
are less dramatic than the change above e ~ 1 GeV/fm^. Therefore, all bound cc states 
could be dissociated at similar energy densities, which lye near 1 GeV/fm^. 

Alternatively, the ip' and the Xc could dissociate above e ~ 1 GeV/fm'^ and the dissociation 
of the J/\l/ could start at e=2.2 GeV/fm^, if we interpret the change in the {J /'^) / DD 
ratio, below and above e=2.2 GeV/fm^, as a step behaviour. In this context, the steep 
drop of the {J/'^)/DD ratio in the bin(s) of largest N (figures |^, ^ |^), cannot be in- 
terpreted in a natural way. The steps of {J/'^)/DD remain to be established through a 
direct measurement of J/\l/ and DD absolute yields as a function of {Ex, N, e). 
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In the above discussed picture, three QGP signatures appear in nuclear coUisions at en- 
ergy density larger than ~ 1 GeV/fm^: 

a) J/\l/ suppression (figure |^ (a)) -which could be due to bound cc states dissociation- 

b) strangeness enhancement (figure ^ (b)), possibly due to equilibration of ss in QGP as 
opposed to hadrons, 

c) the invariant mass m{e~^e~) excess at m below the p mass possibly due to a p 



change and/or to increased production of the lowest mass glueball state in QGP 
This coincidence of QGP signatures, suggests a change of phase at e ~ 1 GeV/fm^ as 
expected |l|]. 

From the above discussion it follows, that a direct measurement of open charm production 
in nuclear collisions appears essential for the physics of the Quark Gluon Plasma phase 
transition. Furthermore, if enhanced over expectations, open charm in nuclear collisions 
defies theoretical understanding. 



5. Possibilities for future measurements 



A measurement of open and closed charm production in Pb+Pb collisions as a function 
of energy below the SPS top energy of y/s = 17 GeV searching for the disappearance of 
the seen J/\l/ suppression in central Pb+Pb collisions at a certain ^/s, could prove clearly 
the QGP phase transition. Using the same nuclei at different ^/s and looking only at 
central collisions, differences due to different nuclear profiles drop out. No currently ex- 
isting or planned experiment at SPS is however able to perform this measurement without 
major upgrades, though one proposal (NA6i) could significantly improve the identifica- 



tion of open charm production through a better determination of the decay vertex 
An upgrade of NA50/NA6i, or alternatively a completely new experiment, could possi- 
bly achieve this goal. The study could also be performed at the Relativistic Heavy Ion 
Collider (RHIC) using lower energy and/or large and small nuclei, and in fixed target 
experiments at RHIC favoured because of higher luminosity as compared to the collider 
mode, important for a low energy scan. 

It would be also important (and easier than the above) to measure the J/"^, DD and DY 
absolute yields per collision, below e=l GeV/fm^, by using the most peripheral (not yet 
investigated) Pb-|-Pb collisions or collisions of lighter nuclei at the highest beam energy 
at SPS (yi = 17,19 GeV). 

An other piece of information important for the understanding of charm production in 
nuclear collisions would be the direct comparison of the {J/'^)/DY and the {J/'^)/DD 
ratios in nuclear collisions at < 19 GeV and mp + p collisions at the Tevatron. Teva- 
tron reaches an energy density similar to or larger than the one estimated in very central 
S+S collisions at 200 A GeV [45|. Therefore it would supply a comparison for these points 



and a continuation of the absorption line fitted through the p+p and p+A data measured 
by NA50 (figure |^), or otherwise. Differences due to the change of dominant production 



13 



mechanisms of charm in pp colhsions as compared to A+B, p+p, can be accounted for 
theoretically. A high Et cut could additionally help in sorting out 'central' p + p colli- 
sions. This comparison should be done possibly in the very same dimuon mass region for 
all processes (also DY) e.g. using Monte Carlo's tuned to p + p Tevatron data. 

This comparison would answer the question, if the energy density is indeed the only criti- 
cal variable for the appearance of a thermalised QGP state with 3 effective flavours u,d,s, 
or whether there is also a critical volume (e.g. as measured by the L variable: V ~ L^). 
Furthermore, at present the comparison of nuclear collisions to p-|-p and p-|-A data is done 
at the same energy and not at the same energy density. This issue is important, since 
if for example the energy density is the only critical scale variable, the QGP should be 
formed also in elementary collisions like pp at a higher beam energy and the same energy 
density. 

Further it is important to search for thresholds in the production of many particles e.g. 
fl, which was found to be enhanced by a factor 15 above p-|-A data in Pb+Pb collisions 
at 158 A GeV |^6[ in the energy density region corresponding to the green stars in figure 
(b). Similarly interesting would be a measurement of the invariant mass of e~^e~ in low 
energy densities. 



6. Conclusions 

In this letter, consequences resulting from the viable possibility that the dimuon invariant 
mass (m(/i"^/i~)) enhancement, measured by the NA50 experiment in the intermediate 
mass region (IMR): m(0) < m(/i"''/i~) < m(J/\E'), in S+U and Pb+Pb collisions at a/s 
19, 17 GeV, reflects a DD enhancement over expectations, are worked out. 

The dependence of the J/\E' and the DD yields per collision in Pb+Pb collisions on the 
mean number of participants has been estimated. This dependence reveals the non ther- 
mal features of charm production at this energy. The ~ A^° '' dependence of the J/\I' yield 
(figure 0) suggests strong dissociation of J/\[' with higher centrality. The dissociation is 
stronger than the absorption seen in any other hadron, e.g. p in Pb+Pb collisions. The 
N dependence of the DD yield of N^'"^ (figure |I|) indicates also non-thermal open charm 
production at this energy, showing up in an excess rather than reduction as compared to 
the thermal expectation. 

If the dimuon excess observed by NA50 is partly or solely due to open charm, it is appro- 
priate to search for an anomalous suppression of J/\E' as compared to the total open charm 
production, rather than to the DY process. We therefore investigated here the {J /'^)/DD 
ratio in Pb+Pb collisions and we find it to decrease approximately as ~ A^^^ (figures |^, ^). 
This is the N dependence expected for the J/\E' if it were completely dissociated in quark 
gluon matter and were later dominantly formed through cc quark coalescence, assuming 
N ~ volume of the cc environment and charm quark multiplicity of one in events with 
charm. In that case, based on coalescence arguments, the [J /^) / DD ratio could be used 
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to estimate the volume of the charm environment, which may reflect partly the size of 
the quark gluon plasma. This is probable under the assumption that the flnal measured 
J/\l/ is dominated by the J/\E' originating from cc pairs which travel through the plasma 
volume, an assumption which may hold only for large plasma volumes, i.e. for the most 
central collisions. 

A further consequence of a possible open charm enhancement is that the J/\l/ over the 
DD ratio appears to be suppressed already in S+U collisions as compared to p+A colli- 
sions, unlike the {J /'^)/ DY ratio (figure The il) / DD ratio would also be additionally 
suppressed as compared to the ^ / DY in both S+U and Pb+Pb collisions. These phe- 
nomena could be interpreted as onset of dissociation of bound charm states above energy 
density e ~ 1 GeV/fm^. 

We estimated and compared the dependence of the {J j^) j DD ratio and of the kaon mul- 
tiplicity per volume (assuming K/N ~ K/V = kaon number density) in several collisions 
and ^/s as a function of the initial energy density. We find that both the kaon number 
density and the ratio {J/'^)/DD exhibit dramatic changes at the energy density of 1 
GeV/fm^, as demonstrated in figure |^ This is the main result of this paper. 

It follows that three major QGP signatures {ss enhancement, p changes and J/\l/ suppres- 
sion) all appear above the energy density of ~ 1 GeV//m^, which is the critical energy 
density for the QGP phase transition according to lattice QCD. 

This discussion underlines the importance of a direct measurement of open charm produc- 
tion in nuclear collisions, and of other experimental investigations proposed in section 5, 
for the understanding of ultrarelativistic nuclear reactions and the dynamics of the Quark 
Gluon Plasma phase transition. 
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Figure 1. DD yield per collision in arbitrary units in Pb+Pb collisions at 158 A GeV, as 
a function of the number of participating nucleons N. The line shows the result of fitting 
the function / = ciV". See text for the fit results. 




Figure 2. J/\l/ yield per collision in arbitrary units in Pb+Pb collisions at 158 A GeV, as 
a function of the number of participating nucleons N. The line shows the result of fitting 
the function / = cN°'. See text for the fit results. 
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Figure 3. Ratio of J/\l/ to DD in arbitrary units in Pb+Pb collisions at 158 A GeV (closed 
points) and in S+U collisions at 200 A GeV (open points), as a function of the number of 
participating nucleons N. The lines show the result of fitting the function / = cN" to the 
S+U (line above) and to the Pb+Pb data points (line below). See text for the fit results. 
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J/ VDDbar=(J/ VDY)/(N°-'') 
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Figure 4. Ratio of J/^/DY to A^o-45±o.ii Pb+Pb collisions at 158 A GeV as a function 
of the number of participating nucleons N. The 9 closed star points correspond to the 
N values at which the DD yield was measured. The lines show the result of fitting the 
function / = cA^" to the closed stars. See text for the fit results. 
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Figure 5. The open points show the ratio of J /'^ to Drell Yan {DY) in p+A, S+U and 
Pb+Pb colhsions normahsed to 200 A GeV, as a function of the path of J/'^ through 
nuclear matter L pO|JT9| . The closed points show the ratio of J/\l/ to DD in S+U and 
Pb+Pb collisions at 200 A GeV estimated here. The open stars show the L dependence 
of the J/'^ /DY and the J/^! /DD ratio in p+A collisions. All J /^ /DD ratios shown in 
this figure are in common arbitrary units. 
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Figure 6. The J /'^ / DD (measured/'expected') ratio is shown as a function of the initial 
energy density (e) achieved in the coUisions investigated. 
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Figure 7. (a) The J/'^/DD (measured/ 'expected') ratio is shown as a function of the 
initial energy density (e) achieved in the colhsions investigated, (b) The kaon (~ K'^) 
multiphcity over the effective volume {V = (n-A- R^^^J ■ (\/l2- -R^ong)) of the particle source 
at thermal freeze out, in the center of mass frame, is shown as a function of the initial 
energy density (e). The above effective volume is smaller than the real source volume but 
proportional to it. 



